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We study the singularity created in the supercritical collapse of a spherical massless scalar field. 
We first model the geometry and the scalar field to be homogeneous, and find a generic solution 
(in terms of a formal series expansion) describing a spacelike singularity which is monotonic, scalar 
polynomial and strong. We confront the predictions of this analytical model with the pointwise 
behavior of fully-nonlinear and inhomogeneous numerical simulations, and find full compliance. We 
also study the phenomenology of the spatial structure of the singularity numerically. At asymptot- 
ically late advanced time the singularity approaches the Schwarzschild singularity, in addition to 
discrete points at finite advanced times, where the singularity is Schwarzschild-like. At other points 
the singularity is different from Schwarzschild due to the nonlinear scalar field. 



PACS number(s): 04.70.Bw, 04.20.Dw 
I. INTRODUCTION 

The spherically-symmetric collapse of a scalar field has 
been studied extensively in the last few years. It was 
shown that when one changes the value of a parameter 
which characterizes regular initial data (while freezing 
the other parameters), there exists a critical value of the 
parameter beyond which an apparent horizon (and con- 
sequently a black hole) forms. For parameters smaller 
than the critical ones, no stable configurations form, and 
the field disperses to infinity jl],[|] . Despite the significant 
advance gained recently in the understanding of dynam- 
ical gravitational collapse, the spacetime interior to the 
apparent horizon has received little attention. In partic- 
ular, the nature of the spacetime singularity predicted 
to form inside the black hole in supercritical collapse, 
remained unresolved. 

The linearized scalar-field perturbations of the inte- 
rior of the Schwarzschild spacetime have been known 
for two decades now ||. This leads naturally to the 
following questions: To what extent is the spacetime 
singularity in nonlinear collapse similar to the linearly 
perturbed Schwarzschild singularity? As the backre- 
action of the scalar field is expected to change the 
Schwarzschild singularity, what are the features and 
properties of this singularity? The goal of this paper is to 
study these questions within the model of the fully non- 
linear collapse of a minimally-coupled, self-gravitating, 
spherically-symmetric, massless, real scalar field. We 
shall construct a simple analytical model to describe the 
singularity which evolves in this model, and confront it 
with fully-nonlinear numerical simulations. 

The organization of this paper is as follows. In Section 
H we present the analytical model and its main conclu- 
sions. We find a generic solution for the pointwise be- 
havior at the singularity. This singulari ty i s spacclikc, 
scalar polynomial, and strong. In Section III we describe 
the numerical simulations and compare these simulations 



with the predictions of the analytical model. We find full 
compliance between the analytic results and the numer- 
ical simulations. We also study numerically the spatial 
structure of the singularity. We make some concluding 



remarks in Section [V 



II. ANALYTICAL MODEL 

The linear perturbation analysis of the interior of the 
Schwarzschild spacetime with a massless scalar field pre- 
dicts the scalar field to diverge like lnr along radial 
t = const lines ||. Here, r is the radial Schwarzschild 
coordinate defined such that circles of radius r have cir- 
cumference 2nr (r is timelike inside the black hole) and t 
is normal to r (t is spacelike inside the black hole). (There 
is a gauge freedom in t — see below.) Therefore, near the 
singularity at r — one would no longer expect the linear 
solution to be valid. Rather, one would expect the back- 
reaction on the metric to be significant. Consequently, in 
order to describe the singularity in this model one would 
need to solve the fully nonlinear Einstein-Klein-Gordon 
equations. Moreover, in the case of dynamical collapse 
the formation of the singularity is a nonlinear effect, and 
no linear perturbation analysis can be made in the first 
place. 

The analytical solution of the full Einstein-Klein- 
Gordon equations turns out to be very difficult even in 
spherical symmetry. However, the analysis of linearized 
perturbations of the Schwarzschild singularity ^ finds 
that at large values of t the interior regions tend to a 
state which depends only on r. For any finite r near the 
singularity the linear spherical scalar field vanishes like 
i~ 3 at large t. One would expect, therefore, the devia- 
tions from Schwarzschild to be small at t — ► oo (at finite 
r). It looks reasonable, therefore, to look for a simplified 
analytical model where dynamical variables depend only 
on r also in the nonlinear case, at least for large t. 
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Let us consider the general homogeneous spherically- 
symmetric line-element 



ds 2 = h(r) dt 2 + f(r) dr 2 + r 2 dVt 2 



(1) 



where dQ 2 = dO 2 + sin 2 9 d(j> 2 is the metric on the unit 
2-sphere. The field equations are given by 

(fr + f 2 -f)/(fr 2 )=<i>> 2 (2) 



(h'r -hf + h)/ (hr 2 ) = $' 2 (3) 



(h 2 )' (f i r /') + 2rf (hh" \h' 2 ^ - 2fh 2 

/ (Arfh 2 ) = -& 2 (4) 



and by the Klein-Gordon equation V AI V M< 1 > = for the 
scalar field <f>. Here, a prime denotes differentiation with 
respect to r. Because of the spherical symmetry and 
the homogeneity, the Klein-Gordon equation is readily 
integrated to 



d>'(r) = d y/\f/(hr*)\, 



(5) 



where d ^ is an integration constant, which will be 
shown below to be arbitrary, and to describe a pure gauge 
mode. We next use Eq. (|J) to eliminate the scalar field 
from Eqs. (§)-(§). We obtain 



h'r hf ■ h ■ i —7 



2h r { fj 2/ 



- 2d 2 4 = 0. 
r 4 



(6) 



(7) 



(8) 



We assume a formal series expansion for /, h, and seek a 
generic solution for a spacelike singularity. We find that 

f(r) = -09 + 1) Cr^ 2 - £±|^±4cV^ 



+ O (r^+ 6 ) 



(9) 



Here, f3 > — 1 and C > are free parameters. This solu- 
tion suggests the the complete formal series expansion is 
of the form 



f(r) 
h(r) 



1 

d 2 /r 2 



E 



4 (AC) 

< (A C) 



„(/3+2)n 



(12) 



It can be shown that the expansion coefficients a^ h can 
be found uniquely for any n [B. For the scalar field we 
find 



$(r) = v^+Tlnr + ^a* (/?, C) 



.(/3+2) 



(13) 



We note that because the Klein-Gordon equation is linear 
in the scalar field whereas the Einstein equatio ns M )- 
([|) are quadratic in the overall sign in Eq. (fLl|) is 
arbitrary. 

We next discuss the genericity of the solution. The no- 
tion of a general solution for a system of nonlinear differ- 
ential equations is not unambiguous. However, one can 
count the number of free parameters in the solution, and 
compare it with the expected number. If these two num- 
bers are equal, than the solution is generic, in the sense 
that it captures a volume of non-zero measure in solution 
space. From the physical point of view one would expect 
two free parameters (one because of Birkhoff's theorem 
and one due to the scalar field). Apparently, our solu- 
tion above has three free parameters, namely, d 2 ,C,f3. 
However, the parameter d 2 reflects a pure gauge mode. 
That is, one has the freedom to re-scale the coordinate t 
by t — > t' = T(t). In such a gauge transformation d 2 is 
changed (and in general the parameter d 2 can be replaced 
by any smooth function of t) without changing the phys- 
ical content of the solution. (The parameter d 2 remains 
a constant parameter which does not depend on t for the 
class of transformations t — > t' = pt + q for constant p 
and q.) Thus, the scalar field ( |l3| ) and curvature scalars 
below are independent of d 2 . Consequently, our solution 
depends on two physical degrees of freedom (/3 and C), 
and is therefore generic. 

The curvature can be described, say, by curvature 
scalars. We find that to the leading order in r 



R 



-Q3+4) 



C 



(14) 



Mr) = d 2 Cr' 3 + d 2 f3{(3+ l } c 2 r 2 P+ 2 + O (r^+ 4 ) (10) 
09 + 2) 2 V 7 



/ VJTT (5/3 2 + 12/3 + 6) o , , 

$(r) = JJ+l lnr + — K — 5— >-Cr p+2 

(/3 + 2) 3 

+ O (r 2 ^+ 4 ) . (11) 



7? RM^ _ 

f 1 ^ — (J2 



„-(2/3+8) 



1 ^iivna lh ^ 



2f3 2 + 2(3 + 3 
(/3 + if C 2 



-(2/3+8) 



(15) 



(16) 



Here, R is the Ricci scalar, R^ v is the Ricci tensor, and 
R^vpa is the Riemann-Christoffel tensor. From these ex- 
pressions the singularity is clearly scalar polynomial. We 
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FIG. 1. The Penrose diagram of the simulated space- 
time. Prior to the characteristic hypersurface the geometry 
is Minkowskian. Due to the supercritical initial pulse a black 
hole is formed. For an outgoing null ray at tio the spacelike 
singularity is hit at v — v*. Future null infinity is at v = oo. 

next show that it is strong in the Tipler sense || , namely, 
any extended physical object will unavoidably by crushed 
to zero volume upon hitting the singularity. Let us de- 
note by x° the coordinate tangent to the worldline of an 
object inf ailing along a radial t = const worldline, and by 
r its proper time, set such that r = at r — 0. Then, the 
(00) tetrad component of the Ricci tensor, to the leading 
order in r, is given by 

*(0)(0)M = 8^^t- 2 . (17) 

By Ref. , Eq. (|l7|) implies that the singularity is strong 
in the Tipler sense, as 

f &t> [ dr"R mo) (r") (18) 

J To J TO 

diverges logarithmically at r = 0. 

Finally, let us define a 'tortoise' coordinate r* by 
9r,r t — —gtt- Then, null coordinates are defined by 
t = (v — u)/2 and r* = (v + it)/2, with ingoing v and 
outgoing u. For a certain outgoing ray u = uq = const 
which runs into the spacelike singularity, the latter is hit 
at some finite value of advanced time v = u* (see Fig. [I]) . 
(Future null infinity is located at v = oo.) We then find, 
to the leading order in — v, that 

r(uo, v) = [d 2 /{P + 1)] 1/4 («* - v) 1 / 2 . (19) 

We note that as this model is homogeneous (the met- 
ric functions and the scalar field depend only on r), it 
actually describes just the pointwise behavior of the sin- 
gularity. In general, the singularity is expected to depend 



not only on r but also on t. This dependence on t cannot 
be captured by our homogeneous model. The inhomo- 
geneities will be studied therefore numerically. 

III. NUMERICAL SIMULATIONS 
A. Initial data 

We constructed a numerical code based on double-null 
coordinates and on free evolution of the dynamical vari- 
ables 0. The code is stable and converges with second 
order (see below). We specify regular initial data on the 
characteristic hypersurface, which is located at it = Uj 
and v — Vi (see Fig. [l]). Prior to the initial pulse the 
geometry is Minkowskian. Then, due to the nonlinear 
scalar field, the geometry is changed, and an apparent 
horizon forms. At late advanced time (v — > oo) the ap- 
parent horizon coincides with the event horizon (EH), 
and the black hole approaches asymptotically the (exter- 
nal) Schwarzschild black hole Q. 

The details of the collapse model, the formulation of 
the characteristic problem, and the description of the 
code are given in Ref. 0. Here, we shall describe them 
very briefly. The null coordinates u, v are defined such 
that v = r and u = (r — ro)/r u o on the character- 
istic hypersurface, where ro = r(ui,Vi). (Note that 
this coordinate v is not identical to the coordinate v 
in Eq. (|l9|).) The initial data we specify are the fol- 
lowing. The scalar field vanishes everywhere on the 
characteristic hypersurface, except for the compact in- 
terval v\ < v < t>2, where the scalar field has the form 
<&(«,-, v) = ^4sin 2 [7r(w — Ui)/(i>2 — t>i)]. The rest of the 
initial data are chosen such that the constraint equations 
are satisfied on the characteristic hypersurface. The pa- 
rameter which we vary is the amplitude A in u). 
For small values of A the initial data are subcritical, and 
no black hole forms. For large values of A, however, the 
initial data are supercritical, and a black hole forms. We 
shall discuss below only the latter case, namely, the inter- 
nal structure of the black hole which forms in supercriti- 
cal collapse. In the simulations described below the free 
parameters in the initial data are chosen to be ro = 10, 
r„ = -1/4, V! = 10, v 2 = 16, and A = 0.29. (These 
initial data correspond to a Minkowskian geometry for 
v < v\ or u < Ui.) These initial data result with a black 
hole with final mass Mr w 4.02. The value of Mf can 
be easily determined from the value of r at the appar- 
ent horizon: r at the apparent horizon is just twice the 
mass of the black hole. Figure || displays the mass of the 
black hole as a function of advanced time v. For v < 11 A 
the domain of integration does not intersect the apparent 
horizon, and therefore the mass of the black hole is not 
shown. (We note that because our code does not include 
the origin, we cannot find the apparent horizon of black 
holes with infinitesimal mass.) We find results similar to 
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FIG. 2. Black hole mass as function of advanced time v. 
The mass is evaluated by means of the value of r at the ap- 
parent horizon. 

those reported below also for other values of supercritical 
initial-data parameters. (We have also studied the case 
of nonlinear perturbations of a pre-existing Schwarzschild 
black hole, by specifying the appropriate initial data. In 
this case we find that the properties of the singularity 
are similar to those found in supercritical collapse. In 
the latter case, however, the dynamical evolution of the 
spacetime geometry brings out the nonlinear aspects in 
a sharper way, and therefore we shall describe here only 
this latter case.) 

The only changes we introduced in the code which is 
described in Ref. j?]] are related to a careful approach to 
r = 0: First, with uniform spacings of the grid in ad- 
vanced time one would hit the singularity within a finite 
integration time and subsequently the code would crash. 
(However, one still needs to preserve the dynamical re- 
finement of the mesh in retarded time.) We resolve this 
difficulty by checking the expected distance (in advanced 
time) to the singularity on the last outgoing null ray of 
the grid for each advanced-time step by means of the 
above analytical results (||) and (0). We then change 
the advanced-time step accordingly, such that crashing 
into the singularity is avoided, and in principle the code 
can run forever. (In practice, because of the floating- 
point arithmetic limits of the machine the approach to 
the singularity is limited. However, one can approach 
r = sufficiently close to infer the asymptotic behavior 
from the simulation.) Second, one would like to study the 
fields approaching different points along the singularity. 
Namely, for different values of u*. This can be done in 
the following way. In Ref. |7| we introduce a procedure 
called "chopping" , which is nothing but excision of the 
domain of integration slightly beyond the apparent hori- 
zon. Here, as we are interested in approaching the space- 
like singularity at r = 0, this procedure is performed only 



up to some finite value of advanced time v c . The later v c , 
the larger , and the closer we arrive to the vertex point 
in the spacetime diagram |l| which is at v* — > oo. This 
way we can arrive at different points along the spacelike 
singularity, and study the spatial dependence of the fields 
at the singularity. 

B. Numerical simulations of the black hole interior 

We perform numerical simulations of the interior of the 
created black hole. Let us first confront the predictions 
of the homogeneous model described in Sec. |l| with the 
pointwise behavior at the singularity as obtained from 
the fully-nonlinear and inhomogeneous numerical simu- 
lations. Then, we shall also study the phenomenology of 
the spatial structure of the singularity numerically. 

1. The pointwise behavior 

We first discuss the behavior of the metric functions 
near the singularity. In double null coordinates the met- 
ric function g uv is expected from Eq. ( |To| ) to behave near 
the singularity like g uv cx r@ to leading order in r, because 
gtt = 2( duj dt){ dvj di)g uv , and dv/ dt = 1 = —du/dt. 
(The gauge freedom in t is manifested just in a multiplica- 
tive factor, and not in the exponent f3 which concerns us 
here.) Figure || (A) displays F = —2g uv along an outgo- 
ing null ray at u — uq, as a function of r. The straight 
line in Fig. ^ (A) indicates a power-law behavior. The 
asymptotic power-law behavior is checked more carefully 
in Fig. U (B), which shows the local power index of F as 
a function of r. This local power index is nothing but the 
exponent fi in Eq. ([h]). As (3 approaches asymptotically 
a constant value, the power-law behavior is verified. 

The behavior of the metric function r(u, v) is shown 
in Fig. H (A), which displays the distance from the sin- 
gularity in terms of advanced time — v, as a function 
of r along the same outgoing null ray as in Fig. || (at 
u = t*o). Again, the straight line indicates a power-law 
behavior, whose local index is displayed in Fig. || (B). 
From Eq. ( |l9| ) this power-law index is expected to ap- 
proach asymptotically a value of 2. This is indeed found 
numerically, within a numerical error of 5 x 10 -3 . 

Another prediction of the homogeneous model re- 
lated to the metric functions is that g rr /guv oc r 2 , to 
the leading order in r. Let us now find the analogue 
of this relation in the inhomogeneous case, and com- 
pare the homogeneous and inhomogeneous cases numer- 
ically. First, in Schwarzschild coordinates V^r V„r 
is readily shown to equal g rr . In double- null coordi- 
nates, /"V^rV^r = 2g uv r tU r jV = — (4/F)r <u r tV , where 
F(u, v) = —2g uv (u, v). Combining these two expressions, 
one obtaines g rr (u,v) — — F / (Ar jU r jV ) . Consequently, 
we find that g rr /g uv = —l/(2r tU r }V ). Figure || shows 
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FIG. 3. The behavior of the metric function g uv along an 
outgoing null ray: (A) F = —2g uv as a function of r. (B) the 
behavior of the local power index of F as a function of r. 
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FIG. 4. (A) v* — v as a function of r. (B) The local power 
index in (A) as a function of r. 
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FIG. 5. The value of b — —l/(2r iU r, v r 2 ) as a function of r. 
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FIG. 6. The scalar field $ as a function of r. 
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FIG. 7. The metric function F as a function of r for three 
values of the grid parameter (N = 20, 40, 80); Dashed-dotted: 
N = 20, dashed: N = 40, and solid line: N = 80. 
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FIG. 8. The scalar field $ as a function of r for three val- 
ues of the grid parameter (N = 20,40,80); Dashed-dotted: 
N = 20, dashed: JV = 40, and solid line: N = 80. 



6 = — l/(2r iU r. t ,r 2 ) along the same outgoing null ray at 
Mo, as a function of r. This ratio approaches asymptoti- 
cally a constant value, which confirms this prediction of 
the homogeneous model. (The gauge freedom in t, which 
is translated here into a gauge freedom in both u and v, 
changes here just the value of this constant.) 

The behavior of the scalar field <I> along u = uo is shown 
in Fig. |[ The straight line here indicates a logarithmic 
divergence of as implied by Eq. (prf). 

The stability and convergence of the code are demon- 
strated in Figs. [?] and [|. Figure [?] shows a magnified 
portion of Fig. ^ and Fig. ^| shows a magnified portion 
of Fig. ^ for different values of the grid parameter N, 
namely, for N = 20,40, and 80, where N is the (initial) 
number of grid points per a unit interval in both the in- 
going and outgoing directions. The behavior in the two 
figures clearly indicates stability and convergence with 
second order. 

The homogeneous model also implies that (3 appears in 
two places, which numerically are independent. Namely, 
in the exponent of h(r) [given by Eq. ([!(])] and in the am- 
plitude of the scalar field [given by Eq. jll])]. Figure |§ 
shows the relation of these two values for (3 for two diffcr- 
rent outgoing null rays: The top panel shows the values 
of j3 for the two null rays as determined from Eq. Jl(i| ) 
as functions of r, and the bottom panel the ratio of the 
two expressions for (3 for the two outgoing rays, again as 
functions of r. We chose here the two outgoing null rays 
to be such that for one of them (3 is positive and for the 
other (3 is negative. The lower panel in Fig. ^ indicates 
that the two expressions for (3 are asymptotically equal 
to within one part in 10 4 . 

The last prediction of the homogeneous model we test 
with the numerical simulations is the behavior of the 
Ricci curvature scalar given in Eq. (114) . The top panel 



in Fig. |l0| displays the behavior of the Ricci curvature 
scalar R along the same outgoing null ray as the data in 
Figs. ||-|| as a function of r. The value of R can most 
easily be evaluated from the right hand side of the Ein- 
stein equations: The Einstein-Klein-Gordon equations 
are R^ = 2$^$^, and the Ricci scalar is simply given 
by R = 2g^ v <& ,^3>,i/, which in double-null coordinates and 
spherical symmetry reduces to 



R = 

r 



(20) 



In order to evaluate the expression ( |l4{ ) for R we need 
first to express C . This can be done by comparing our 
result that g rr — —F/{Ar u r^ v ) with the leading order in 
r expression for / from (^|). We find that to the leading 
order in r 



C 



F 



1 



-08+2) 



4/3+1 r u r „ 



(21) 



The bottom panel of Fig. shows the ratio of R as cal- 
culated from Eq. ( pp[ ) and R as calculated from Eqs. ( |l4| ) 
and @, as a function of r. This ratio asymptotically 
deviates from unity by 3 x 10~ 4 , despite the growth of R 
by a factor of order 10 20 . 



2. The inhomogeneous behavior 

The inhomogeneous behavior near the singularity can 
be demonstrated by the variation of (3 along the singu- 
larity. Namely, (3 is no longer a parameter, but (3 = f3{t). 
[We would also expect C = C(t).] Another issue which 
can be studied numerically is the overall sign issue of the 
scalar field. This sign is of importance, as its meaning is 
related to the divergence of $ to +oo (for a minus sign) 
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FIG. 9. Top panel: The values of (3 for two outgoing null 
rays as functions of r. Bottom panel: The ratio of (3 as deter- 
mined from the exponent in g uv and as determined from the 
amplitude of the scalar field $ for the same two outgoing null 
rays. 

or to — oo (for a plus sign). We find that there are regions 
where $ diverges to +oo and regions where <E> diverges 
to — oo. Figure [ll] (A) shows (3 as a function of v». Fol- 
low values of (3 is positive. Then, /3 changes its sign, 
arrives at a minimum at which its value is —1, increases 
again, and finally, at late v*, approaches asymptotically 
the value of —1. Although the numerical simulations in- 
dicate that (3 has at least two more local minima, be- 
cause of the numerical noise it is hard to locate them 
accurately. This can be done, however, in terms of the 
scalar field. The sign of the scalar field can be deter- 
mined from c = r<& ri which, to the leading order in r is 
-fVP + 1- Clearly, whenever c(v*) vanishes with a non- 
vanishing spatial gradient (namely, with dc/ dv* ^ 0), (3 
has a local minimum with a value of —1. From Figs. [ll] 
(B), (C), and (D) this happens at v* = 13.4, 22.5, and 
66.4. We note that these values of are not the same 
as the values of v along the EH where the quasinormal 
modes have their extrema or vanish. However, they may 
be related to the quasinormal ringing, as their frequen- 
cies are similar. It is hard to quantify this and establish 
the hypothetical relation between the two oscillations, as 
we do not have enough cycles of the oscillations of c on 
the singularity. Yet, Fig. [ll] suggests that c oscillates a 
number of times before it approaches a vanishing value 
for — > oo. In Fig. [ll] (B), (C), and (D) the abscis- 
sae were chosen to overlap, such that one could infer the 
overall behavior. 

We find, that for a number of discrete points at finite 
values of «*, (3 = —1. From Eq. ( |TT| ) the scalar field 
<fr vanishes at these points. One could therefore ask, to 
what extent the singularity at these points is similar to 
the Schwarzschild singularity. Let us assume, that in the 
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FIG. 10. Top panel: The behavior of the Ricci scalar R 
as a function of r. Bottom panel: The ratio of R as deter- 
mined from the fully nonlinear and inhomogeneous numerical 
simulations and the value of R as determined from the homo- 
geneous model as a function of r. 

fully inhomogeneous case the leading order for the met- 
ric functions and the scalar field are given by the same 
form as the leading terms in the homogeneous case, ex- 
cept that the free parameters are replaced by functions 
of t. (Note that higher order terms will have a different 
form in the inhomogeneous case.) Then, in order that 
the singularity will be Schwarzschild-like at these points, 
both the scalar field and its gradient need to vanish there. 
The gradient of <E> will vanish only if 8(3/ dt vanishes. To 
see if indeed d(3 / dt = at the local minima of (3 let us 
consider the neighborhood of the first minimum of /?(«*)• 
Figure [l2| shows /3 as a function of v r near this first min- 
imum. Figure [l^ suggests that indeed 8(3/ dt = at the 
minimum, and consequently the singularity is expected 
to be Schwarzschild-like there, in the sense that at least 
the leading order for the scalar field and its gradient van- 
ish there, and (3 = — 1 as in Schwarzschild. Of course, 
one would need to consider also higher-order terms in or- 
der to study their contribution to the scalar field and to 
the geometry in an inhomogeneous model for a more de- 
tailed comparison of the singularity at these points and 
the Schwarzschild singularity. 

Another point to be made here is related to the limit 
v* — * oo. At this limit (3 approaches a value of —1, and 
the spatial gradient of (3 vanishes. This is evident from 
the behavior of c for large u* . Consequently, the singular- 
ity asymptotically approaches the Schwarzschild singu- 
larity at late v*. This conclusion agrees wih the conclu- 
sion drawn from the linear analysis in Ref. However, 
for any finite v* (except for the discrete points for which 
the scalar field and its gradient vanish on the singular- 
ity) (3^—1. This means that the scalar field does not 
vanish on the singularity, and by its backreaction on the 
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FIG. 11. (A) P as a function of v, along the spacelike sin- 
gularity. (B), (C), and (D): c = r<I> ir as a function of v* along 
the spacelike singularity, for different sections of the latter. 



geometry the latter is not Schwarzschild. The deviation 
of the geometry from Schwarzschild is expected, however, 
to be small for small (3 + 1, and decrease with growing t. 
In this sense the nonlinear numerical simulations confirm 
(at least in this model) the dictum of Doroshkevich and 
Novikov ||] that "a black hole has hair neither outside 
nor inside." 

IV. CONCLUDING REMARKS 

We remarked above that in order to study the structure 
of the singularity in this model without the simplifying 
assumption of homogeneity, one could perhaps replace 
the free parameters in the solution with functions of t. If 
this were indeed the case, the leading terms are expected 
to retain their form. However, the higher-order terms in 
the solution will not. This is due to the appearance of 
non-vanishing spatial derivatives in the Einstein-Klein- 
Gordon equations. The solution of the full inhomoge- 
neous equations looks like a formidable challenge even in 
the special case of spherical symmetry and a scalar field. 

Another point concerning the scalar field is related to 
its serving as a model for a dynamical physical field. The 
scalar field was introduced as a toy model, because it has 
a radiative mode in spherical symmetry. In the more 
realistic case of gravitational waves the field does not 
have radiative modes in spherical symmetry, and there- 
fore has no non-trivial dynamics. One could ask the fol- 
lowing. To what extent do the spherical symmetry and 
the scalar field capture the essential dynamics of realistic 
black holes with vacuum perturbations, or the dynami- 
cal collapse of spinning fields? Realistic spinning black 
holes are known to have a Cauchy horizon, which is ex- 
pected to transform into a null weak singularity || . We 
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FIG. 12. f3 as a function of u». The circles are the data 
from the numerical simulations, and the solid line is a linear 
interpolation of these data. 

have shown that in the model of the collapse of a spher- 
ical scalar field a spacelike singularity evolves. Will this 
be the case also in the collapse of gravitational waves? 
Although as yet there is no answer to this question, we 
could consider the model of the nonlinear perturbations 
of a spherical scalar field of a Reissner-Nordstrom black 
hole. In this model one also has a Cauchy horizon in the 
unperturbed spacetime, and one can show that indeed a 
spacelike singularity evolves PJlf],^, [TT| , p^ | . It turns out 
that there are some important differences between the 
spacelike singularity we find here (with vanishing elec- 
tric charge) and the spacelike singularity in the model 
with charge, despite the similarity of the formal series 
expansions. First, the parameter /3, which we find here 
to be greater than —1, is restricted to positive values 
in the presence of charge. Second, we find in the un- 
charged case that in the limit v* — * oo the parameter 
/3 — ► — 1. In the charged case (3 grows very rapidly for 
large i;*, and it is reasonable to expect it to diverge as 
Vif — ► oo. This difference in the asymptotic behavior of 
f3 can perhaps be attributed to the following. In the 
uncharged case as u* — > oo the outgoing Eddington-like 
coordinate u e — > — oo. However, because of the pres- 
ence of the Cauchy horizon in the charged case, the limit 
V* — > oo corresponds to a finite value of u e . The infinite 
range of retarded time in the uncharged case acts to in- 
finitely redshift the fields, an effect which does not have 
an analogue in the charged case. 

Finally, one might worry about the role that the scalar 
field plays in this model. Scalar fields are notorious for 
being related to peculiar phenomena. For example, scalar 
fields are known to destroy the oscillatory nature of the 
BKL singularity (r|. Therefore, even though we find the 
singularity with a spherical scalar field (both in the un- 
charged and charged cases) to be monotonic, one perhaps 
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could not infer from that on whether the (as yet hypo- 
thetical) spacelike singularity inside spinning black holes 
is monotonic or oscillatory. 
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